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Abstract. - Current statistics can be calculated in various ways. Event-based approaches use 
the statistics of the number of events occuring during a given time. Time-based approaches use 
the statistics of the time needed to reach a given number of events. By analyzing normal as well 
as anomalous statistics of nonequilibrium currents through a two level system in contact with two 
different reservoirs, we investigate the conditions under which these different statistics do or do not 
yield identical predictions. We rely on the continuous time random walk formulation introduced 
in our earlier work [Phys. Rev. E 77, 051119 (2008)]. 
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Introduction. — The study of currents is of fun- 
damental importance in statistical mechanics. Average 
energy and particle currents are central in determining 
whether a system is in equilibrium. Indeed, at equilib- 
rium the net probability current between any pair of sys- 
tem states is zero, i.e., detailed balance between all pairs 
of states is satisfied. This of course means that all net 
currents through the system are zero. External driving 
mechanisms bring a system out of equilibrium and may 
lead to non-zero currents. Examples of driving mecha- 
nisms include contacts with reservoirs of different temper- 
atures or chemical potentials, and external forces. The 
resulting departure from detailed balance is responsible 
for an associated non-zero entropy production [T]. 

The study of current statistics of systems out of equilib- 
rium has become an important tool to analyze small sys- 
tems where fluctuations around average behavior can be 
significant. In macroscopic systems only average currents 
are of interest because fluctuations are too insignificant to 
be observable. However, in small systems the fluctuations 
of currents around the average provide rich information 
about the nonequilibrium dynamics of these systems, as 
evidenced by the considerable recent interest in the field 
of photon and electron counting statistics [2]-[6]. Further- 
more, current fluctuations have recently been connected 



to the generalized version of the Second Law of thermody- 
namics via a variety of fluctuation theorems [7HT2] . How- 
ever, there is a fundamental aspect of current statistics 
that has, to our knowledge, not yet been fully explored, 
namely, the circumstances under which ergodic conditions 
are or are not satisfied. Ergodicity involves the equivalence 
between ensemble-averaged and time-averaged statistical 
properties at equilibrium. Here we extend the notion of 
ergodicity to nonequilibrium currents by exploring con- 
ditions under which event-based and time-based current 
statistics are asymptotically equivalent. 

It is of particular interest to address these questions in 
the presence of anomalous statistics [13] that are known 
to significantly affect the ergodic properties of the sys- 
tem [MHTT] . Effects of anomalous statistics were first 
observed in averaged current measurements on macro- 
scopic disordered semiconductors and amorphous solids, 
and were described in terms of continuous time random 
walks (CTRW) [T8H20| . Nowadays, anomalous statistics 
can be measured by photon counting experiments on sin- 
gle nanodevices such as blinking single quantum dots |21j , 
single nanocrystals [I4j, and single molecules [22] ■ Since 
these experiments can give access to the full probability 
distribution of the fluctuations it becomes appropriate to 
extend the concept of ergodicity beyond the average to 
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higher moments. 

Definitions and main results. To describe the 
issue, we introduce a system in which certain elementary 
events occur as time t progresses, for instance, the net 
transfer of a mass or of a charge in or out of a system. 
The process is stochastic, that is, there is a distribution 
of the number of events occurring in a given time interval, 
or there is a distribution of times at which one observes a 
given number of events. The current associated with such 
events is usually defined as I = k/t under the assumption 
that the number of events, fc, scales linearly with time. 
We will call this the "normal" scaling, and generalize the 
definition of the current to recognize the possibility that 
the number of events does not grow linearly with time. 
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Normal behavior thus corresponds to a = 1, and we 
call other scaling behaviors "anomalous". In particular, 
we focus on anomalous behavior with sublinear scaling, 
< a < 1. As we shall see, this scaling is consistent 
with the CTRW formalism and, more specifically, with the 
subordination principle of so-called anomalous statistical 
processes [5DH5^ . 

The scaling information is not sufficient to describe the 
full statistics of the current. We need to be more pre- 
cise about the nature of the distribution that defines the 
behavior of the current and current moments. The tradi- 
tional approach to current statistics is to count the number 
of events that occur during a given time interval. Mathe- 
matically, this means that one treats A: as a random vari- 
able and i as a parameter. The associated probability 
distribution is Pt(fc), the probability that k events occur 
in time t. The current moments according to this point of 
view are calculated as 



(2) 



where m is a positive integer. The average current is given 
by the first moment, {I)t- The subscript serves as a re- 
minder that time here is a parameter. 

An alternative approach is to consider time to be the 
random variable and the number of events the parameter. 
The associated probability density is Pk{t), the probability 
that a time t is needed to observe k events. The current 
moments according to this point of view are calculated as 



^ dt t°"^ Vk{t)' 
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With this approach, we use a different symbol for the cur- 
rent (J) simply to stress the difference. The average cur- 
rent is now given by the first moment {J)k- We stress that 
([2]) and (|3|) require different types of measurement of the 
current. The first measures the current during a fixed time 
interval t and the second measures it until a given num- 
ber of events k has occurred. The statistics is of course 



lim ^— ~ lim , 



obtained by repeating the experiment many times or by 
doing it simultaneously on independent copies of the sys- 
tem. 

At this point it might be tempting to say that the cur- 
rent statistics are ergodic when asymptotically the statis- 
tics of events and the statistics of times lead to identical 
results, 

' k 

that is, 

lim = lim (J")fc. (5) 

To our knowledge this concept of ergodicity would be new 
for two reasons. First, it defines an equivalence of time 
and ensemble averages at the level of the currents instead 
of ordinary system observables. The time average of a 
current cannot be expressed in terms of the fraction of time 
that the system spends in a given state along a trajectory, 
as do conventional approaches to ergodicity (e.g., [15]). 
Furthermore, currents are not uniquely specified by the 
system states and require that the reservoirs responsible 
for the transitions be specified as well. Second, we have 
not restricted this statement to averages {m =1) but have 
stated it for all moments (m > 1). 

One of our two main results of this paper is the follow- 
ing: 

The ergodic condition (jlj is satisfied for normal 
statistics, a = 1, but fails for anomalous statistics, 
< a < 1. 

Our second goal in this work is to introduce current 
statistics for which an ergodic condition is satisfied for the 
entire range < a < 1. This is achieved by dealing with 
the current even more directly, so that time as a random 
variable governed by the probability density Vkit) enters 
through the inverse moments, 
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An alternative ergodic condition would then be the asser- 
tion that asymptotically the statistics of events and the 
statistics of inverse times lead to identical results, 



lim A = lim fc^^t-"™)*,, 



(7) 



that is. 



lim - lim (r")fe. (8) 

t— >oo fe— >oo 

The second main result of this paper is the following: 

The ergodic condition ([7]) is satisfied for all values 
< a < 1. 

While we could describe and justify these assertions in 
a very general way for currents that scale as ([T}, it is 
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more instructive to do so using a simple model that high- 
lights the principal issues. Our model system is sketched 
in Fig. [1] and falls within the framework considered in 
Ref. 

The model. — The system consists of two levels, i — 
and i = I, connected to two "particle" reservoirs that we 
call left (L) and right (R). The nature of the "particles" 
does not matter other than that the current is associated 
with the flow of these particles between L and R through 
the two-level system. The nature of L and of R also does 
not matter other than that L and R must of course be 
different in some way. It is the asymmetry between the 
two that generates the nonequilibrium constraint. We 
have chosen the simplest possible system for our example, 
namely, one that allows the presence of either no parti- 
cles or one particle but no more than one particle. Such 
a model can be viewed as a single-site asymmetric sim- 
ple exclusion process [24] and has been used, for example, 
to model electron transport in a single-level quantum dot 




Fig. 1: Two-level system embedded between a left (L) and a 
right (R) particle reservoir. 

Transitions between the states of the system occur be- 
cause of particle flow in or out of the system to one or 
the other reservoir. We follow the CTRW formulation 
in [23] (also called a semi-Markov process, e.g., [261127]) 
to describe these transitions. The transitions between the 
states are triggered by one or the other of the reservoirs, 
and since these are distinct and correspond to different 
mechanisms, we label the transitions with an index v — 
L, R. Suppose that the system arrives at state i' at a given 
time and that its next transition is to state i at a time t 
later via mechanism v. We assume the distribution 7/;^^, (t) 
for this transition to be separable. 



(9) 



Here Pqi and P^q are the transition probabilities down- 
ward and upward, respectively, triggered by reservoir v, 
and the waiting time distributions ipoit) for transitions 
from i = to « = 1 and ipi (t) for transitions from i = 1 to 
i = are assumed to be independent of v. 

To explore the consequences of anomalous statistics on 
various notions of ergodicity, we choose waiting time dis- 
tributions that decay asymptotically as 



The times tq and ti are characteristic times (but not in 
general first moments, which in fact diverge when a < 1). 
As we shall see shortly, this scaling is consistent with the 
scaling ([1]) of the current. The case a = 1 is representative 
of waiting time distributions that decay at least as fast as 
and these have a first moment r^. In Laplace space, 
to lowest order in s the distributions behave as 



Ms) = 1 - {nsT- 



(11) 



We note that a direct connection to thermodynamics can 
be made for a = 1 by, for instance, identifying the two 
reservoirs as heat baths of equal temperatures T but dif- 
ferent chemical potentials fi^, and imposing the conditions 
iPwn)/iPoiTo) = expifi^/ksT) 123]. A detailed discus- 
sion surrounding our choice of waiting time distribution 
can be found in [33] . 

Current statistics. — Having described our setup, 
we are now ready to calculate the various statistical quan- 
tities needed to prove our statements about ergodicity. 
First, consider event-based statistics. We define the num- 
ber of events to be the difference between the number 
of — > 1 and 1 — > transitions (denoted A^oi and A'^io) 
triggered by the v reservoir by time t, i.e.. 



(12) 



It is always understood that a transition occurred at time 
t = 0. In Ref. [33] we used a generating function formalism 
to obtain an analytic expression for the asymptotic form 
of the moments of Vk^ (t). Those results immediately lead 
to the result for the left hand sides of Eq. (j?]) [or ([S])] and 
Eq. © [or ©], 
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F(to + 1) 
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where 



r>iy pi/ 



(13) 
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To calculate the right hand sides, we next consider the 
time-based statistics, that is, the statistics of the times 
needed to reach a given number of events. The probabil- 
ity Pk^ (t) is the probability that k^, events have occurred 
by time t and is simply a time convolution of A^io factors 
ipoit) and A^oi factors ?/'i(i) in some order. The order is in- 
tertwined with the events triggered by the other reservoir, 
but does not matter when we go to Laplace space. Recall- 
ing that convolutions in the time domain lead to products 
in Laplace space, we get for the Laplace transform of P^^ 
the product 



PkAs) = [Ms)f^''[M'S)]''°' 



(15) 



Mt) - (iM)(Vn)" 



(10) 



The numbers of events A'^io and A^oi are random variables. 
However, this can be simplified at long times because 
whereas the total number of transitions 2N = iVgi -I- Niq 
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is large, the difference is small, jiVoi — A^iol < 1, the con- 
straint coining from the fact that no more than one par- 
ticle is allowed in the system. Therefore at long times we 
can set Nqi « Niq k. N, so that 

pK{s) = [Ms)Ms)f- (16) 

In the long time (small s) limit, Eq. (|lip then leads to 



which can be approximated by 

PkAs) « exp[iVln(l-(ro"-HTf)s") 

( , 
= exp s 



(17) 



(18) 



In the last step we have set k^, « N{Pqi — Pfg) and have 
used Eq. (fTi)l . 

Equation (fT8|) is the inverse Laplace transform of the 
one-sided Levy distribution. 



l/a 



r 



l/a -| 
t 



(19) 



The fractional moments of the Levy distribution for < 
a < 1 and q < ct are [281129] 



y du Caiu) 



r(i - g/g) 
r(i-'?) ' 



from which it follows that 

r(i-g/«) 



r(i-9) 



qjo 



(20) 
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The higher positive moments (g > a) are singular. These 
moments allow us to calculate the right hand side of of 
Eq. ^ [or Eq. ([S])]. In particular, the first immediate 
conclusion is that for < a < 1 the moments in the de- 
nominator of the right hand side of Eq. ^ diverge. This 
means that the right hand side of this presumptive ergod- 
icity condition vanishes, whereas the left hand side does 
not, cf. Eq. ([131). Thus, the condition (gl [or indeed 
fails for anomalous statistics. 

On the other hand, when a = 1 one finds that for any 
value of q 

<">'--(^)'- <^^' 

and hence we find for the right hand side of Eq. Q that 



lim , , 



(23) 



This is identical to the result obtained from Eq. (|13p when 
a = 1. Thus, this condition is satisfied for normal statis- 
tics. 



Finally, we turn to the right hand side of Eq. Q or ([51). 
The negative moments can be obtained from Eqs. (|20p and 
([2Tt with q <Q. We immediately get 



r(TO + l) a™ 



(24) 



Using this result, we find that the time-based current 
statistics leads to the moments 

Vim 4- 1 ~) 

hm (C)fc„ = hm C(i""")fe. = 
fc„->oo k^,^oo i (am -f- Ij 

The last expression is identical to that occurring in 
Eq. (IT^ . This then proves the ergodic condition ^ or 
(0) for all < a < 1. 

Conclusions. — In this letter we have considered the 
current statistics of systems away from equilibrium in 
which an external constraint produces a current. The 
waiting time distribution between successive elementary 
processes that give rise to the current decay as t~°'~^ at 
long times, with < a < 1. The scaling is "normal" when 
a — \ and "anomalous" when < a < 1. The elemen- 
tary events may, for example, be a transfer of particles or 
of charge. Specifically, we have addressed the question of 
the equivalence between different ways of calculating the 
current statistics. 

Two setups can be used to directly measure a (properly 
scaled) fluctuating current / = k/t" . One can fix the time 
interval t and measure the number k of events occuring 
during this time or, alternatively, one can fix the number 
of events k and measure the time t required for this num- 
ber of events to occur for the first time. We have shown 
that in both of these setups, the average current as well 
as the higher moments of the current are asymptotically 
the same, that is, limt_>oo(fc'")/r'" = lim/o^oo 
for ah < a < 1, cf. Eqs. (0), and A simi- 
lar conclusion, but for a different model involving average 
particle velocities in tilted periodic potentials, has been 
reached in [34l . When the current moments are calculated 
indirectly using the average time statistics cf. 
Eq. ([31), difficulties arise in the case of anomalous scaling. 
This occurs because this approach is in fact an incorrect 
way of calculating current moments except for normal scal- 
ing, where all three calculations of the current moments 
yield the same asymptotic results. 
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